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Abstract 

We explicitly show that split-universal extra dimension (split-UED), a recently suggested exten- 
sion of universal extra dimension (UED) model, can nicely explain recent anomalies in cosmic-ray 
positrons and electrons observed by PAMELA and ATIC/PPB-BETS. Kaluza-Klein (KK) dark 
matters mainly annihilate into leptons because the hadronic branching fraction is highly suppressed 
by large KK quark masses and the antiproton flux agrees very well with the observation where no 
excess is found . The flux of cosmic gamma-rays from pion decay is also highly suppressed and 
hardly detected in low energy region (Ey < 20 GeV). Collider signatures of colored KK particles 
at the LHC, especially qıqı production, are studied in detail. Due to the large split in masses of 
KK quarks and other particles, hard pr jets and missing Er are generated, which make it possible 


to suppress the standard model background and discover the signals. 
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I. INTRODUCTION 


Recent observations by PAMELA |1| and ATIC/PPB-BETS bh consistently suggest a 
new primary source of energetic cosmic electrons and positrons in the neighborhood of our 
solar system, since high energy electrons/positrons lose their energy quickly within about 1 
kpc. These resutls have attracted enormous attentions and stimulated many interpretations 
of the primary electric source from both astrophysics and particle physics, for example, 
the nearby pulsars and dark matter decay i or annihilation |6]. The existing data of 
cosmic electrons and positrons cannot distinguish among these possibilities, and we expect 
that future data from cosmic-ray and accelerator experiments such as Fermi and Large 
Hadron Collider (LHC) will give us some clues on which interpretation is more promising 
and ultimately correct. 

In this paper, we consider one of the most attractive models of dark matter, universal 
extra dimension (UED) with its minimal realization (mUED)|7Z,|8] and a recently suggested 
variety split-UED |9]. When the reflection symmetry (dubbed Kaluza-Klein parity) about 
the mid point of the orbifold extra dimension is exact, the lightest Kaluza-Klein particle 
(LKP), which is odd under the parity operation, is absolutely stable and is a good candidate 
of dark matter [10]. It turns out that in most cases of universal extra dimension scenarios the 
LKP is the first Kaluza-Klein (KK) photon (Bı), and a pair of charged leptons is produced 
in the annihilation process of Bı. Interestingly enough, the mass of dark matter suggested 
by the ATIC/PPB-BETS data (~ 600 — 700 GeV) exactly coincides with the one that gives 
the right_relic density for the LKP dark matter after taking co-annihilation channels into 
account {11}. 

On the other hand, the measured data of the cosmic ray antiproton-to-proton flux ratio 
between 1 and 100 GeV by PAMELA follow the expectations from secondary production 
and strongly constrain contributions from dark matter particle annihilation gi Therefore 
it is required to naturally reduce antiproton flux in any dark matter model. Split-UED is 
promising in the sense that production of antiproton is quite suppressed. The cross section 


for dark matter annihilates into a quark pair is 
2 2 252 
(0B, Bı —qqU) x ms, / (Mp T Min) 


where mp, and m,, are the masses of Bı and the first KK quark (qı), respectively, and it 


can be significantly suppressed by increasing the masses of KK quarks, i.e. splitting the 


201708.00323v1 


chinaXiv 


KK quarks from other particles. However we should keep in mind that predictions of the 
observed antiprotons from dark matter annihilation can vary quite a lot for different adopted 
diffusion models as we will see later. As hadronic production is suppressed in split-UED, we 
would also expect much less production of cosmic gamma-ray whose main source is decay of 
pions. This prediction will be tested by forthcoming data from the Fermi experiment 


The split spectrum of Kaluza-Klein particles is realized by introducing a bulk mass term 
for quarks in split-UED keeping the KK parity exact 1 One immediate consequence 
is that the 5D wave functions of quark fields are quasi-localized at the boundaries and 
induce violation of the translational symmetry along the extra dimension. Accordingly the 
KK number conservation is violated in the quark sector, and KK-even gauge bosons could 
directly couple to the zero mode quarks at the tree level, which is forbidden in mUED. 
KK-even bosons could be copiously produced at colliders and it is promising to discover 
them, for example, by the 2nd KK gluon resonance at the LHC |9]. Another consequence 
is the heavy KK quark decay to the SM quark and KK boson, which will generate harder 
jets with large missing momentum compared to that in mUED due to a large mass splitting 
between KK quarks and other KK particles. Such differences may give us some handle to 
tell split-UED apart from mUED at the LHC. 

This paper is organized as follows. In Section [L| we review split-UED and show how 
the spectrum and the interactions are modified compared with those in mUED. Note that 
the renormalization group (RG) running is taken into account in our spectrum calculation. 
We then calculate the annihilation cross section to leptons and quarks with varying “bulk 
mass parameter”, u, which is the only new extra parameter in split-UED. In Section 
we estimate the signals of cosmic-ray_positrons, electrons, and antiprotons from LKP dark 
matter annihilation in both mUED and split-UED, and compare our predictions with 
recent experimental data. Moreover, the diffusive gamma-ray is studied as well. In Section 
IV| we consider the productions of KK quarks and gluons at the LHC, and we also discuss 


the possibility of discovering these signals and mass reconstructions. We finally conclude 


our studies in Section 


1 In warped geometry KK parity could be kept in a different way hd. 
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II. MASSES AND COUPLINGS IN SPLIT-UED 


Our set-up is based on the all SM fields are “universally” propagating in one extra dimen- 
sion, which is compactified on an orbifold, St/Z2, with the boundary points y = +L (the 
orbifold radius R = 2L/7). The extension we have introduced is an scalar background ®(y) 
that couples to each colored Dirac fermions V,(x", y) in the 5D bulk where the corresponding 
standard model quark field q (=SU(2) singlets u$, dẹ and doublet Qz) resides on the chiral 
zero mode after orbifold projection. A common Yukawa coupling A(= Ae, = Aue, = Aag) 
is assumed for all quarks so that any further source of flavor violation beyond the CKM 
matrix is automatically avoided and the model becomes more predictive ?. Here we neglect 
the family indices but one can easily extend the model to include more generations. Finallly 


the 5D action is given by the form 


L 
Ssplit-UED = fef dy [Cmuen = AB (y) Y, (z, y)V,(z, y)| , (1) 
=} 


where the summation runs for all quark fields q, and Lmurn contains the usual mUED terms. 

If we choose an odd background profile ®(—y) = —®(y) under the inversion symmetry 
about the middle point y = 0, we can see that the KK parity, which transforms the fields as 
(x,y) —> P(x, —y), V(x, y) — £75V,(x, —y), is still a good symmetry of the Lagrangian. 
For quantitative concreteness we consider the simplest mass profile which respects the re- 
quirement: A(®(y)) = we(y), where u is the bulk mass parameter and e(y) is a step function 
defined as +1 for 0 < y < L and —1 for —L < y < 0. As (A — 0) or equivalently (u — 0), 
split-UED is continuously reduced to mUED thus we call this limit as “mUED limit” below. 
For gauge bosons, scalars, and leptons (l), their KK decompositions are the same as those 
in MUED, while for quarks, their odd masses will affect the KK decompositions. The de- 
tails of the all KK decompositions in Eq. (I) are presented in Appendix A. Notice that in 
the case of KK decomposition, we also label the KK parity in addition to the KK number 
n = {n7 (n odd), n*(n even)}. 


The tree level mass for the n-th KK fermion is given by 


mire = \/p? + k2, (2) 


m tree = n/R, (3) 


In 


2 In general the Yukawa coupling can be assigned for any fermion field in the bulk and they can be different 
from each other. 
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if we neglect their masses coming from electroweak symmetry breaking. Here kn is deter- 
mined by the boundary conditions which are given by: 
Kn- 


— |u| tan kn L, (4) 


We choose the minimal universal boundary conditions at the energy cutoff scale A that all 
boundary kinetic terms vanish °, and i expression for the one-loop correction to 


the n-th KK fermion masse is given by |1 


9 G 
1+ = (= C3 (Nac) log AR 


where ag = g&/4m and we neglect the Yukawa corrections in this formulat. The sum 


= tree 
M fa — M Fa 


, (6) 


is over all SM gauge groups under which the fermion f is charged, and CẸ(N) is the 
quadratic Casimir operator of the fundamental representation N in gauge group G, e.g., 
G MN = (N? — 1)/2N and we take C IN = Y? with Yp being the hypercharge of 
f. If we choose (a1, a2, a3) & (0.010, 0.033, 0.094) at Q? = (620GeV)?, colored particles will 
get bigger corrections (~ 14% for up, dh and 16% for Q} ) than leptons (~ 1.9% for l} and 
3.1% for IŁ). 

In Fig. [I] we plot the spectra of 1st KK quarks and leptons including 1-loop corrections. 
When pu = 0, i.e. the mUED limit, quarks are quite degenerate with leptons even we take 
the QCD corrections to the masses into account. However when we increase the bulk mass 
parameter u, the KK quark masses become larger and larger so that we get the split spectra 
as we request. Note that small deviations in quark spectra come from different SU(2) and 
U(1) charge of ur, dr and Qz. In our study, these differences can be neglected. 

The couplings between KK gauge bosons and fermions determine the most interesting 
phenomenological features of the model. All 4D effective couplings could be calculated by 
integrating out the 5D profiles of the relevant fields along the extra dimension. Since 5D 
profiles are the same for different KK gauge bosons (quarks) with the same KK number, the 


, over the SM one is the 


ratio between a general gauge boson-quark-quark coupling 9.4, —~qn—q/ 


3 One can also consider the non-universal boundary conditions as an extension. See for example, Ref. lg. 

4 For the KK top quark, if we consider the large top Yukawa coupling, positive contribution to the tree level 
mass term will cancel the negative one from the radiative corrections, and the overall mass corrections 
will not be affected too much. 
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FIG. 1: (upper plot) The spectra of 1st KK quarks and leptons in the unit 1/R with the given 

bulk mass parameter u. The lepton masses are independent of u in our scenario. Here we included 

l-loop correction for all the masses. The most important correction is QCD correction and it 

is roughly 14 ~ 16% for quarks. For leptons the correction is small (~ 2%). (lower plot) The 

magnified figure of the upper one in wR = (3.60,3.90). SU(2) doublets (tz, Qz) get the larger 


1-loop correction than singlet quarks (ur, dpr) get. 


same for fixed KK numbers, where m, n and l are KK number of A, q and q’, respectively. 
We introduce a function Gmni(wR) to parameterize how such a ratio depends on the bulk 


mass parameter u times the radius R: 
JAm—qn—qı/ 90 = Gran uR) , (7) 
where go is the SM gauge interaction between Ap, qo and qj. Among those couplings, the 
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FIG. 2: The coupling constant ratio g4,—q,—-q/9o = Giio vesus wR. The dependence comes from 
the change of the 5D profiles. Notice that the coupling ratio first increases to its maximal value 


around 1.14 for wR ~ 1 and then decreases monotonically. 


most interesting one in phenomenology is the coupling between the 1st KK gauge boson 
Aj, the 1st KK quark qı, and SM quark go, which is crucial on LKP pair annihilation and 
productions of qı and the 1st KK gluon gı at the LHC. The dimensionless function G110(u R) 
is calculated from the 5D profiles of Ai, qı and qo given by Eq. BD, 68), G9), and (41) in 
Appendix A, and is plotted in Fig. 2las a function of wR. At u = 0 limit we get the mUED 
result, which is the same as the SM coupling. As u gets larger, the 5D profiles are localized 
so that the resultant overlap among the fields get smaller. 

Taking care of all the corrections and effects of u we can now calculate the branching 
fractions that LKP pair annihilates into quarks and leptons. At the non-relativistic limit 
(8 = ,/1—4m%,/s — 0), the cross section for LKP pair annihilates into fermions can be 


well approximated as: 


2 


291Cr + Ov), (8) 


= 2 
9mm Bı 


GuoluR) 
(1+ ri 


where a useful numerical term Cp = NY, + Yi) is introduced with the color factor Ne = 
3(1) for quarks(leptons). The mass ratio between the LKP and KK fermion is parameterized 


by a parameter rf defined as r; = mj /Mp, and we ignore the small mass corrections coming 


T 
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FIG. 3: The ratio of LKP pair annihilation into quarks and charged leptons vesus uR. The mUED 
limit corresponds to u = 0 and the ratio of cross section to quarks to leptons is about 47% but it 
quickly goes down below 10% as u > 1.3/R. The shaded region (o4q/oll > 10%) is not preferred 
by PAMELA antiproton data. 


from SU(2),, U(1)y gauge group so that my, = Mf ~% Mfp: In MUED limit, the ratio 
between productions of the lepton and the quark is mainly controlled by the ratio of Cy from 
the hypercharges and color factor of fermions. For leptons and quarks the corresponding 
terms are `, C, ~ 3.18 and $` q Ca = 1.90, respectively. Obviously the mUED does have a 
sizable hadronic annihilation channel and a carefully study of antiproton flux is important. 
In the case of split-UED, one can easily adjust r, to suppress the hadronic annihilation cross 
section. In Fig. B] we plot the ratio of annihilation cross section into quarks and charged 
leptons ogq/o,7 with wR in the range of (0,5). We find the ratio goes down very quickly as 
the bulk mass parameter goes larger since o ~ 1/(1 + re). By turning on the bulk mass 
parameter we can easily avoid the antiproton excess in the cosmic ray detection experiments. 
We also note that the soft gamma-ray from decay of hadrons, especially from pions, is under 


control in this set-up. 
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TABLE I: The branching fraction of the LKP pair annihilates into various final states, for different 
qı mass, and different bulk mass parameter u. The numbers shown are summed over generations. 
We neglect the W-boson final state because the SU(2) component of LKP is too small. Here we 
fix the LKP mass to be 620 GeV and include the radiative corrections to all the masses according 


to Ref. |17]. 


H (GeV) 


0 200 400 600 800 


BR(BıBı > qq) 


BR(B,B, > I) 
BR(B,B, > v?) 
BR(B,B, > ¢¢*) 3.1% 


II. COSMIC POSITRON, ANTIPROTON AND PHOTON 


LKP dark matter Bı will mainly annihilate into a fermion pair while the W-boson and 
Higgs boson modes are suppressed, as shown in Table[l] For the quark-pair final state, due to 
the QCD hadronization process, a bunch of hadrons will be produced and sequentially decay 
into positrons/electrons, protons/antiprotons, photons and neutrinos. Moreover, positrons 
and electrons also originate in the leptonic channels: muon and tau (also generate photons) 
decay and direct production in B,;B,; — e*e~ process. In our calculation, we use PYTHIA 

8| for the simulation of the QCD hadronization and the decays of particles. Since these 
stable particles from the dark matter annihilation in our Galaxy will propagate to our solar 
system and be observed as cosmic-ray signals, we estimate the predictions from our model, 
the split-UED, and compared with experimental data in this section. Due to the fact that 
excesses were observed in the positron fraction in PAMELA [1] and the flux of electron plus 
positron in ATIC/PPB-BETS ldla while the ratio of antiproton to proton, p/p, is consistent 
with the astrophysical background, we first explain positron and electron data from the dark 
matter annihilation. Then base on the parameters set by fitting the ATIC/PPB-BETS and 
PAMELA data, we calculate the p/p and photon flux. Since the W-boson and the Higgs 
boson productions are highly suppressed in UED models we are considering, we will not 


take into account their contributions in the following calculations. 
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Ko(kpce?/Myr)|L(kpe) 


TABLE II: The parameters in the diffusion models which are comparible with B/C and generate 


minimum (M2), median (MED) and maximal (M1) positron flux. 


A. Positron and Electron 


In this section, we calculate the cosmic-ray positrons and electrons from the BB, annihi- 
lation which account for the PAMELA and ATIC/PPB-BETS data. Since the dark matter 
annihilation generates the same amount of electrons and positrons, we only show the formula 
for positron flux below. After being produced from the dark matter annihilation process, the 
positrons will propagate in the magnetic field of the Milky Way. Because the magnetic fields 
are tangled, the motion of the positron can be described by a diffusion equation. Neglecting 


the convection and annihilation in the disk, the steady state solution must satisfy 
~ , 0 P 
K(E) V’ fa l(E, T) + a BE) (ET) + QE, F) = 0, (9) 


where fe+ is the number density of e+ per unit kinetic energy, K(E) = Ko(E/GeV)? is the 
diffusion coefficient, b(£) = 10716 (E/GeV)? is the rate of energy loss and Q(E,7) is the 


source of producing e*. In the Bı dark matter annihilation case, 


Q(B, #) = 3 (42) Dew, (=P=), (10) 


where dN,+/dE is the energy spectrum of et obtained by using PYTHIA, the index 7 runs 
over all quark and charged lepton pairs and p(7) is the dark matter profile. In our numerical 
calculations, we adopt an overall boost factor Br and the isothermal halo model [19| which 
is given as ° 


\ Po 
Phalo(7) = T+ (r/r.?’ (11) 


where r = |r] is the distance from our Galactic center, re = 3.5 kpc and po is the parameter 


that is adjusted to yield a dark matter local halo density of 0.3 GeV/cm? [19] in our solar 


5 Effectively we are using p?(r) = Brp?.,),(r) 


10 
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system. Considering the diffusion zone of a cylinder with radius R and half-height L, the 


solution for the positron flux at the Earth can be written in a useful form 


c 
D2M(E) = fa (E, ro), (12) 
where c is the speed of light, 
margt as a ; Caro \n. /mr MT 3 Gao 
fet (Ero) = f pr tk 2 Qnm E \Jo ( R sın (=) exp (oF) -+ S X 
(13) 


where EB is the energy in a unit of GeV, Jo is the zeroth-order of the first kind Bessel 
function, ro ~ 8.5 kpc is the distance from Milky Way center to the Sun, Cn» is the n-th zero 


of the function Jo, 


x= et (14) 
and : E 
2 Cr . (MZ = 
Qnm(E) = sone | arrJ( SE) f desin( TEE- 2) QEP (18) 


The Jı in the Eq. (l5) is the first-order of the Bessel function and the parameters ô, Ko are 
set in the simulations of diffusion models to coincide with the observed cosmic-ray data, 
especially the Boron to Carbon ratio (B/C) . The values for different diffusion models 
we adopted are listed in table [I] |22). 


In addition to e* fluxes from dark matter annihilation, there exist secondary e* fluxes 
from interactions between cosmic-rays and nuclei in the interstellar medium. We use the 


approximations of the e~ and et background fluxes [23 |24 


rim 0.16711 i 5 j i 
oe (BE) = TPO EED GeV~*cm~ “sec” tsr >, 
0.7697 
se E) — = eV 1 2 1 1 
ow) 1+110F15 + 600E? + 580 E42 ee 
sec 4.5 E07 i 2 1 1 
er (E) GeV “cm “sec” “sr, (16) 


~ 14 650E23 + 1500542 


where F is in units of GeV. Therefore, the fraction of positron flux and total flux of positron 


plus electron are 
ODM + Of 


re 17 
PPM + PDM 4 pse + ROE + pse u7) 


Lt 


201708.00323v1 


chinaXiv 


1000 C T T T T T T TT T T T T T T TIT T T 
F m, = 620 GeV, boost factor = 200 
1 


L propagation model: MED 


È 100 4 
s S = = 
g background my 
= J 
“lu ATIC J 
i Torii et al 2008. ey || 


DuVernois et al. 2001 
Torii et al. 2001 

— mUED 

— split-UED (KK quarks decoupled) 


°¢ > @ nu 


quark mode 


ji L L L L i | L L 
10 100 1000 
E (GeV) 


FIG. 4: The flux of electron plus positron from LKP dark matter annihilation using MED diffusion 
model, compared with experimental data pi pi bd pi Note that only lepton mode is inclulded in 
the split-UED, since all the KK quarks are decoupled. 


and 


PPM + BOM + oie 4 kT + Hee, (18) 
respectively, where k is a free parameter which is used to fit the data when no primary 
source of et flux exists lod La, For our numerical simulation, we take the mass of dark 
matter, mp, to be 620 GeV. In Fig[4] we show the total flux of electron and positron with 
the background from Eq.(i6) and recent experimental data Ld p bd bd. In order to have 
a better fit for the data, a boost factor of 200 is needed, which is consistent with the value 


(200) chosen in_Ref.|2}. Boost factor is known to have origins such as local clumps in dark 


matter ale al also ladh, Sommerfeld enhancement effect by a long range attractive 


more origins. In our case, without assuming a new attractive force or further tuning of mass 


force |3 , BA, the Breit-Wigner type resonance effect and possibly many 


spectra to get the large resonance, we tend to assume that boost factor mainly comes from 


clumps. We also show in Fig[4] the contributions of quark mode and leptonic mode from 


the dark matter annihilation. The e* from quark mode are much softer than that in the 


leptonic mode. This is because that the main source of e* in quark mode is the hadron 


12 
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FIG. 5: The same as Fig[4] but only for split-UED, and different diffusion models, M1, MED and 


M2, are shown together. 


cascade decay while in leptonic mode there exists a direct production of B,B, — ete, and 


e~ from u™ and t* decay are harder. We can easily see that the peak is mainly contributed 


by the leptonic mode and sharp drop-off is due to the direct production of e*, therefore, 
the predictions of mUED and split-UED are quite similar to each other in the peak region 
since the difference between these two models is just in quark sector. We also compare 
different diffusion models in Fig [5] for split-UED case. MED and M1 models produce more 
soft electrons and positrons than the M2 model does, thus the distributions of the former 
two models are flater. For energetic positrons and electrons, all of three models are quite 
the same with each other. 

For positron fraction, we show in Fig[6] the predictions of mUED and split-UED, using 
MED diffusion model, with PAMELA data Al Since there are more soft et in mUED from 
the quark modes, the distribution is falter compared to that of split-UED. If we adopt M2 
and M1 models, as shown in Fig[/] the fittings seem worse compared the MED case for split- 
UED especially in the energy region around 10 GeV. However, we have decoupled all of the 
higher KK quarks in this plot, i.e. hadronic final states are turned off in the dark matter 


annihilation. With finite mass of a first KK quark, more soft positrons will be produced, 
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FIG. 6: The positron flux fraction predicted from the LKP dark matter annihilation using MED 


diffusion model, compared with PAMELA data [1]. 


therefore, the M2 model can also be consistent with the data. For M1 model, we expect 
that by tuning the parameters, e.g. boost factor or the normalization of the background, it 
would be possible to explain the PAMELA data as well. So, we conclude here that the dark 
matter in both mUED and spilt-UED can be the equivalently good source of the excesses 
in PAMELA and ATIC/PPB-BETS data. 


B. Antiproton 


The propagation of antiprotons through the Galaxy is similar to that of positrons. The 


diffusion equation of antiprotons can be described as 


K, 7° folT, T) — 2 [(sign of 2)Ve(z) fa(T, T)] — 2hô(2) ann fp(T, r) + Q(T, 7) =0, (19) 


where fp(T,r) is the number density of antiproton per unit energy, T is the kinetic energy 
of antiproton and K, = Kop3(p/GeV)° is the diffusion parameter. V.(z) in the second term 
of Eq.(£9) is related to the convective wind that tends to push antiprotons away from the 
Galactic plane, and is assumed to be a constant . Again, like the case in positron, the values 


of Kop, ô and V; for different diffusion models are set to agree with the observed cosmic-ray 
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FIG. 7: Same as Fig[6] but only for split-UED, and different diffusion models, M1, MED and M2, 


are shown together. 


data, especailly the B/C, and their values are listed in table Bi The third term of 
Eq. (19) represents the annihilation of p with the interstellar proton in the Galactic plane, 
h is the half-height of plane which is set to be 0.1 kpc in our calculations. The antiproton 
annihilation rate Pann is given as [36 

Dann = (ny + 4 nue) (Top (20) 


where ny ~ 1/cm? and nye ~ 0.07ny are the number densities of hydrogen and helium, 


respectively; og,” (T) is the annihilation cross section and vz is the velocity of antiproton. 
The form of opp” (T) is given by [3 
ae 661(1 + 0.0115T~°"4 — 0.948799!) mbarn , for T < 15.5 GeV; 
r= (21) 


S67 ->> mbarn for T > 15.5 GeV. 


The solution of the interstellar flux of antiproton in the vicinity of solar system is 38 
Up = 
S(T) = E f(T, ro) (22) 


=e 4 | Mr 
— i >an (2 E 
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Kop(epe? /Myr)| (kpo) |Ve(km/s) 


TABLE III: The parameters in the diffusion models which are comparible with B/C and generate 


minimum (MIN), median (MED) and maximal (MAX) antiproton flux. 


where the dN;/dT is the spectrum of antiproton which we simulate with PYTHIA, and the 


index 7 runs over all quark final states in the Bı dark matter annihilation processes, and 


~ T: VoL Yn(L) 
G= Pa Jo (GR) exp Ea A, sinh(S,,L/2)’ (23) 
with 
4p t v-a]. 
Yal L) = zor | ardor] R) | dzexp ar sinh(S,(L — z)/2)p(r)", (24) 
and 


An = 2hPann + Ve + Kp(T)Sn coth(S,L/2), pa 
V2 40? 

_— c on 26 

2 KAT) FR an 


However, we have to take solar modulation into account to estimate the flux of antiproton 
obtained at the the Earth, which is important for the low energy antiproton. The flux of 
the antiprotons at the top of Earth’s atmosphere can, therefore, be written as ad 


2 
2Mpl TOA i a T oA 


pIvA(T = 


®* (Tis), (27) 


where Troa = Tig — ġr with pr being the solar modular parameter which we take to be 500 
MV in our calculation. In order to compare with the observed data of antiproton to proton 
ratio, we need to include the astrophysical background. Here we adopt simple fittings for 


background antiprotons oy and protons prs provided in Ref.|41 , 


oy = 2 SU [GeV "m?s tsr7}], (28) 
P [4+ 300-1 + 0.0823 


0.9t-1 


———— E 
8 + 1.1¢71-85 + 0.8 ¢1-68 


b ee ae eee | 
a [GeV m?s "sr". (29) 
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The predicted ratios of antiproton to "id ld aw MUED and split-UED are shown in 
Here we take a common boost factor 


Fig [8] with experimental data Lid la Lid [ed 
(Br = 200) which we have obtained from _ PPB-BETS fit. Since the hard electronic 
signals are mainly from local sources within a few kpcs but antiprotons can come from 
distant sources, one may take different values of boost factor for positron and antiproton in 
one’s calculation. Here, partly for predictability and simplicity and partly in order to follow 
common assumption that we live at a typical place in our galaxy, we take a common boost 
factor for all particles in our calculations. The mUED agrees well in the low energy region 
E X 10 GeV, but starts deviating significantly from the recent PAMELA data when energy 
of antiproton becomes higher, if MED diffusion model is adopted. Therefore, mUED seems 
to be disfavored by PAMELA data in the high energy region, and it agrees with PAMELA 
observations only when MIN diffusion model is chosen. However, one can easily see that 
how the situation can be improved in the split-UED case. By enhancing the mass of qı to 
be twice of the l4, the prediction of p/p is reduced by more than a factor of two, and will be 
much suppressed if qı becomes heavier, as shown by the blue and magenta lines in Fig[8] We 
also note that if we adopt MIN model, the split-UED case agrees with PAMELA very well 
even for mg, = 2m. The lower bound of the mass of qı in split-UED can be estimated by 
the upcoming data of PAMELA in the higher energy region, since a small bump is predicted 
at E ~ 200 GeV, and the height of deviation from background is controlled by the mass of 


qı- 


C. Gamma-ray 


For diffusive gamma-ray, there are galactic and extragalactic contributions from the an- 
nihilation of LKP dark matter Bı. The flux of the gamma-ray from the extragalactic origin 


is estimated as |46 


dJ, BrE? p up dN. (1 + z)78/2¢7#/2max 
2 yY Bı C y 
le Z| = otis <a>: in (5 1 z) i 
e 4nmp, HQ i ((1 +2)E) Jä (1 + z)-3 
M 


(30) 
where Qp,, Qm and Qa are the density parameters of Bı, matter (including both baryons 


and dark matter) and the cosmological constant, respectively; pe is the critical density; Ho 
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FIG. 8: Antiproton to proton ratio predicted in mUED and split-UED models, together with 
experimental data. In mUED, both MED (red line) and MIN (green line) diffusion models are 
shown. For split-UED, the cases of mg, = 2m), (blue line) and m,, = 3m, (magenta line) are 


presented, using MED diffusion model. 


is the Hubble parameter at the present time; z is the redshift, and z,, = mg,/E — 1 since 
the maximum energy of photon in the rest frame of annihilation is EF = mp,; the term 
eW7/mex with zmax(E) ~ ae takes the optical depth into account Lag) For 


the numerical results, we use |4 
Qp, h? = 0.1099, Nuh? = 0.1326, Na =0.742, p, = 1.0537 x 10 5GeV/cm?. (31) 


On the other hand, the gamma-ray flux from the annihilation of Bı in the Milky Way 


halo is 
dJy E? 1 dN, 


BE? SNS Laa 2 AVIA 32 
| a) An 2m}, dE (fo W) 7 (32) 
where Phalo iS the density profile of dark matter in the Milky Way, ( des P (Üa) is the 


average of the integration along the line of sight (los). 
Since there are still big uncertainties and ambiguities for modeling the Galaxy center, 
we integrate over the whole sky except for the zone of the Galactic plane (i.e. exclude the 


region with the galactic latitudes |b| < 10°). For the background, we use a power-law form 


18 


201708.00323v1 


chinaXiv 


F Extragalactic ~'0, ‘ 
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FIG. 9: The diffusive gamma-ray predicted in mUED and split-UED models, together with 
EGRET data. The solid lines are for signals with background, blue line is for mUED; green, 
magenta and red lines are for mg, = 2M}, Mg = 3M), and mg, — ov, respectively. Signals of dark 
matter annihilation are shown separately in dashed lies for leptonic final state (red), hadronic final 


states in mUED (blue) and in split-UED with mg, = 2m), (green) and Mmg, = 3m ;, (magenta). 


adopted in Ref. Lat 


d 
es ~ 5.18 x 10 7E7049 GeVem “sr ‘sec, (33) 
bg 


where EF is in units of GeV. 

The total flux of diffusive gamma-ray is shown in Fig] and the signals from the hadronic 
and leptonic final states in B,B, annihilation are presented as well. We should emphasize 
that we do not try to fit EGRET data lid in the plot, since Fermi is expected 
to have more precise results soon, however, we still show EGRET data for reference. We 
notice that the extragalactic contribution from cosmological distance is very small even we 
enhance it by a factor of 10° to estimate the effect of subhalo [46]. The gamma-ray in the 
decay of 7° from hadron cascade decay is softer than that from the decay of T, therefore 
the distribution become flater when more quarks are produced in Bı annihilation. By 
suppressing the fraction of the quark mode in the final state, i.e. comparing mUED (blue 


lines) and split-UED (green, magenta and red lines), the starting point of deviation from 
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the background will shift to high energy ê and the predicted peak at energy about 200 GeV 
will become smaller. The bump at the high energy is a generic prediction for split-UED 
even all of the higher order KK quarks are decoupled, because gamma-rays still come from 


the 7 decay, and it can be further examined by the Fermi experiment in the near future. 


IV. COLLIDER SIGNATURE OF SPLIT-UED 


Having light colored particles below a few TeV in split-UED, the LHC can produce lots 
of KK quarks and gluons via QCD interactions. As it is already shown in Sec. [| one of the 
features of our model is that KK quarks are split from the other particles, therefore, will lead 
to collider phenomenology quite different from the mUED models. For example, the lst KK 
quark, where we take m,, ~ 2mp, through calculations, mainly decays to the 1st KK gluon 
gı and a SM quark, generating at least one high pr jet due to the big mass gap between 
qı and gı. Another one is the existance of tree-level KK number non-conserved interactions 
between KK-even gauge bosons and the SM fermions, as a result, the production cross 
sections of 2nd KK gauge bosons can be substantial. It has been shown that the signals of 
a 2nd KK gluon production can be dominant over the SM background in dijet events after 
imposing a certain set of cuts at the LHC [9]. In this section, we will focus on the scenarios 
of 1st KK quarks and gluon, which can only be produced in pairs. 

The production cross sections involving qı at the LHC depend on the bulk mass parameter 
u. Among the production cross sections of ag”, ag and nag, the o(qiqi) is dominant 
one as the valence u and d partons can contribute to the initial states. These productions 
proceed mainly through the t-channel diagrams with the 1st KK gauge bosons exchanged. 
The most relevant couplings for these productions are qı-q-Vi, where V; = gı, B4, the 1st 
KK W-boson W; and the 1st KK Z—boson Z,. The cross section o(u u,), o(uid,) and 
o(d,d,) are plotted in Fig. [0] with varying u, and we neglect the ng” productions, since 
they are insignificant. For small u, the production cross sections increase with increasing 
u due to the enhancement from the couplings, as shown in FigB]| since the cross section is 


proportional to (Gi19)*. For large jz, on the other hand, those cross sections decrease very 


ë The preliminary result of Fermi shows the data is consistent with the estimated background up to 
E ~ 10 GeV, so by enhancing the mass of qı in split-UED seems to be prefered, if the preliminary result 
is thoughtful. 
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FIG. 10: The cross sections and the mass of wz, as functions of the bulk mass u with fixing 


1/R = 620 GeV. Here, uı(dı) includes both uz; and ur; (dz; and drı) contributions. 


quickly with increasing u (increasing m,,) as expected. For the gıqı and gig; productions, 
the cross sections, which are not shown, are at the same order as that of q,q,. Here all cross 
sections are calculated using CalcHEP by modifying the mUED model file in Ref. : 

The signature of split-UED qıqı productions followed by qı — giq — BiqX is two 
high pr jets plus soft jets/leptons and large missing momentum. The kinematics and the 
decay branching ratios are very similar to those of the corresponding supersymmetric model 
(SUSY) processes, namely the g¢ productions followed by g > gq — X?qX with squark mass 
mg = Mq, gluino mass mg = Mg, and the lightest neutralino mass myo = mg,. Therefore, 
we mimic split-UED collider signatures of qı pair productions at the LHC by scaling the total 
cross sections of SUSY signatures with the model point whose mass spectrum is similar to 
that of the split-UED. Furthermore, because collider signatures of qıqı are the same as that 
of gq, they share the same SM background. Since the SM background for squark and gluino 
productions at the LHC is very well studied , we apply the same cuts in our calculation. 
The model points we study are listed in Table |LV| in which we use ISAJET for SUSY, 
and we take 1/R and u to be 620 GeV and 700 GeV, respectively, for split-UED so that 
Mg © 2mp,. The events are generated by using HERWIG kd for SUSY processes, then we 


regard them as the split-UED events, and detector simulations are carried out by AcerDET 
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1347 GeV|lūz ,dz|1355, 1358 GeV 


1304 GeV 
1263 GeV 
799 GeV 


622 GeV 


TABLE IV: Mass spectrum of split-UED at 1/R = 620 GeV and u = 700 GeV. The corresponding 
SUSY point generated using ISAJET is also listed. 


o after standard cut|| Meg > 1 TeV| Meg > 1.5 TeV 


TABLE V: The acceptances for the split-UED q1q1, gig, and gig; productions after imposing basic 
event selection cuts and harder Meg cuts for the model point in Table Note that these 


numbers are obtained from SUSY events. 
b. Table [V] shows the acceptances of the events after imposing the following basic cuts: 
Meg > 500 GeV, ETmiss > max(100 GeV, 0.2Mer), 100 > 1, nso > 4, (34) 


Meg is the scalar sum of the pr of the first four leading jets; ETmiss is the transverse missing 
momentum; n100 (nso) is the number of jets with pr > 100 (50) GeV. Forthermore, we also 
use harder cuts on Meg with Meg > 1, 1.5TeV. 

We can see that the acceptance is about 40 % for qı pair production after the basic cuts. 
On the other hand, the efficiency to select the gı pair production is much smaller, especially 
when a large Meg is required. For example, the number of events from the Ist KK gluon 
pair production reduces by factor of 1/30 after imposing the basic cuts with Meg > 1 TeV. 
This is because the probability of having high pr jets and large E-pmiss is very low for the gı 
production events at our point. Thus, we expect gg, production is not promising due to the 


very small efficiency found in the simulation. The g;q; production may be easier to detect 
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FIG. 11: Meg distributions for “g1g1”, “q1gı” and “qıqı” from left to right obtained by correspond- 
ing SUSY channels gg, gg and qq. We neglect the y-axis here since we focus on the shapes, and we 


see that the distributions shift to higher energy region from left to right. 


compared to the gg, case with the help of the high pr jet from qı — qg1, although the signal 
and the background separation would be worse compared with qıqı production. Therefore 
we only consider qıqı productions in the following studies. For completeness, we show the 
Meg distributions for 9191, ¢191 and qıqı in Fig. Again these results are obtained from 
SUSY events mimicking split-UED with the same kinematics. 

The total production cross section of q,q1, including uyj,dz1, Uri, dpi, is 7.64 pb, and we 
expect 7640 events for 1 fb~!. According to Table [I] we expect 7640 x 0.37 = 2830 events 
left under the standard cut with Meg > 1TeV. We should also consider lepton veto before 
using the background studied in st Although it depends on the lepton branching ratio, 
we expect more than a half number of events, i.e. 1400, remain with the lepton veto. The 
number of SM background events with the same cut for 1fb~! is less than 300 according to 
the Mee distribution shown in Ref. [54]. Therefore the q,q, signal distribution is well above 
the background distribution for Meg > 1000 GeV. 

Given the enough statistics, we now discuss the possibility to determine some of the 
split-UED particle masses. Fig. shows the Mr» distribution of the two highest pr 
jets under the standard cut with Mog > 1000 GeV. The Mrz is usually used to determine 
the masses of unknown particles which are pair-produced and decay identically. In our case, 
we have two qı’s being produced, and both decay into gig. Here the Mrz is calculated 


from the two highest pr jets, jı and jg and a missing transverse momentum defined as 
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FIG. 12: The Mrz distribution for qıqı production events under the standard cut with Meg > 1 
2 
TeV for 1 fb~+. We can see that the end point of Mr distribution is given by Mmg, — “a ~ 880 GeV. 


May 


The bin size is 10 GeV/bin. Note that we do not apply any lepton veto cut for this plot. 
PTmiss = PTj — PTj and the formula is 


Mr = min [max{ Mr (¢ri, Pj, Miria), Mr(¢re, Pjo, Merial) }] , (35) 


PTmiss=9T1+9T2 


where gr; and gro are two dummy parameters that make up Prmiss, and the minimization is 
taken for all possible sets of qrı and qr2; pj,(j2) is the momentum of j1(2); Mtrial is the trial 
mass that represents the unknown mass of the daughter particle (gı in our process) and Mr 
is the transverse mass. The Mrz depends on the trial mass Miria, and we took Miria = 0 
in Figf2] For our simulation, the two highest pr jets comes from qı — gig, and the end 
point of Mrz should be m,, when the trial particle mass mia) is taken as the true mass of 
gı. Another important fact is that when there are no initial state radiations, the formula of 


the Mrz endpoint for Miria = 0 has a general form given as 


Mest = m, — —, (36) 
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where m4 is the mass of mother particle which is pair produced and mx is the mass of the 
unknown daughter particle, and in our case m4 = Ma and mx = Mg. We can see that 
the Mr» distribution in Fig. [2] has a clear end point which is consistent with the predicted 
value Mo? ~ 880 GeV for the model point in Table[[V]. Although we are unable to measure 
precise masses of qı and gj, we still have a useful information on their combination. 

Finally, let us comment on the comparison with mSUGRA. If we assume the signature 
comes from squark production, with a typical mass spectrum mg ~ 1 TeV and gives the 
similar Mrz end point, the cross section is much smaller than that of the split-UED point. 
Thus, SUSY and split-UED should be distinguished based on the event rates, although the 
kinematical nature of the signal is similar. We can get a rough understanding on the enhance- 
ment of the cross section in split-UED by comparing the helicity structure of squarks/gluinos 
and corresponding KK quarks/gluons 


T See also for similar phenomenon in Little Higgs model. 
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V. CONCLUSION 


The split-UED model proposed in Ref. id is studied in detail. From the basic set-up 
of the model, we calculate all the masses (including one-loop radiative corrections to the 
masses) and relevant couplings for the first KK particles. Characteristic split spectrum of 
KK quark states is realized and enables us to control the hadronic branching fraction of the 
LKP annihilation by one parameter, bulk mass scale (u). In the non-relativistiv limit, such 
a branching fraction decreases very fast by increasing p. 

In order to explain the anomalies from the cosmic ray, we fixed the size of extra dimension 
(1/R ~ 620 GeV) from the peak position in flux of electrons plus positrons observed by ATIC 
and PPB-BETS, with boost factor B = 200. Then all data from PAMELA, ATIC, and PPB- 
BETS can simultaneously be fitted very well. It is important to notice that such a dark 
matter mass in split-UED generically can explain the right amount of relic density of dark 
matter (Q ~ 0.23). To avoid all the strong constraints known up to date, we take mg, Z 2/R 
so that the hadronic branching fraction of LKP annihilation is significantly smaller than the 
leptonic one (043/0 < 10%) and the cosmic antiproton flux is suppressed more than the 
safe rate. The flux of cosmic gamma-rays from pion decay is also highly suppressed and 
hardly detected in low energy region (E, < 20 GeV). However future coverage of high 
energy domain in O(100) GeV will make it possible to detect gamma-rays originated by tau 
lepton which is still sizable. 

At the LHC, the lst KK colored particles (KK quarks and gluons) in split-UED are 
copiously produced and further decay into hard jets, soft jets/leptons and dark matter 
(missing energy). Because of the large mass splitting between the Ist KK quark qı and 1st 
KK gluon gı, we can separate out the production of gig1, gig, and qıqı by using the Meg 
cuts. We focus on the q,q, production,in particular, which has two hard jets so that we can 
distinguish the signals from the SM background. The split-UED q pair signal is simulated 
by scaling the SUSY signatures with the same mass spectrum. After the appropriate cuts, 
our signal is well above the SM background. We calculate the Mrz distribution of the 
qı pair signals, and its end point reflects the information for combination of qı and gı 
masses. Although we are unable to determine the mass of individual particle in split-UED, 
its predictions of the signatures of qıqı production at the LHC, e.g. Meg and Myo, can be 


examined whether they are consistent with the cosmic-ray signatures in the near future, and 
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vice versa. 
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Appendix A: KK decomposition in split-UED 


For all the gauge bosons, scalars, and leptons, their KK decompositions are the same as 
those in mUED. For bosons and fermions that are in the same chirality of the zero mode, 


which has (+, +) boundary condition at |—L, L], the KK decomposition is: 
Bey) = oon(a) ro E sn Fu] oe) +005] | oe) 6T 


While for fermions that are in the opposite chirality of the zero mode, which has the (—, —) 
boundary condition at [—L, L], the KK decomposition is: 


(y) = zÈ cos l $n- (x) + sin Sp] bnt (2) . (38) 


The label n~ or n* here stands for the n-th KK modes with the even/odd KK parity®. 
For the quark, we consider the case in which a SM quark is embeded into Y, = P,W 
component of a 5D Dirac fermion V, where P} = (1+ 5)/2 is the positive chirality projection 
operator. The other case could be considered quite similarly by hd. u by —p. 
14 


convenient way of expressing the KK decomposition is in the form |9,, 


Tiley) = SY  gnt(ll)Xn+(t) + €(Y) Gn (Iyl)Xn-(2), 
W_(x,y) = XO ey) fally nle) + fa- llul) Yn- (2). (39) 


nt n 


The 5D profiles satisfy the following coupled, first-order equations of motion 


0, 


Oy In F UJn = Mn fn 
Oy fn = fin + ™MnGn = 0, (40) 


8 In our convention, the KK number n~ = 2n~ — 1 and n* = 2n* for the {n,n} in Ref. id. 
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with each 5D profile gn+, Gn-, fn+, fn- satisfying the (+, +), (—, +), (—, —), (+, —) boundary 
condition at y = 0 and L, respectively. Once the solution for y C [0, L] is obtained, the 
solution to the whole space is determined from Eq. (89) thanks to the symmetry. 


For the zero modes, Mm, = 0, the equations are separable and the solution is given as 
y 
go (9) = Avexp( | A(0(s)}ds) = Ao exp(uo), (41) 
0 


with the normalization factor Aj = \/p/e?@#4-)). There is no zero mode company for go+ 
but each of the KK modes has its couple and fill the spinor states of the Dirac spinor. The 


solution for the even KK modes is 


Tr ka 
Int = -| + cos(kn+y) + 
L limnt Mnt 


iPS -sintio (42) 


where the KK mass my+ = y/u? + k2, and k,+ = na/L (n € 2). 
To avoid the very light 1st KK quark, we choose u > 0 and the zero mode is quasi-localized 
at the boundary y = +L. For odd KK modes, the solution is 


Li kns 
fn- F | cos(kn-y) a 
m 


ne nT 


sin(kn+y)|, 


sin(kn-y)|, (43) 
where the KK mass Mm, = 4/ u? + k? and kn is the n-th solution of the equation 
kn- = —ptan(k,-L) . (44) 


When u increases from 0 to +00, k,- increases from (n — 1/2)r/L to na/L. In this case, in 


the limit of 44 > +00, all KK modes could be decoupled. 
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